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Abstract 

We present extensive Scanning Tunneling Spectroscopy (STM/S) measurements at low temperatures in the multi- 
band superconductor MgE>2. We find a similar behavior in single crystalline samples and in single grains, which 
clearly shows the partial superconducting density of states of both the -k and a bands of this material. The su- 
perconducting gaps corresponding to both bands are not single valued. Instead, we find a distribution of super- 
conducting gaps centered around 1.9mV and 7.5mV, corresponding respectively to each set of bands. Interband 
scattering effects, leading to a single gap structure at 4mV and a smaller critical temperature can be observed in 
some locations on the surface. S-S junctions formed by pieces of MgB2 attached to the tip clearly show the sub- 
harmonic gap structure associated with this type of junctions. We discuss future developments and possible new 
effects associated with the multiband nature of superconductivity in this compound. 
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1. Introduction 

The surprising discovery by Nagamatsu et al. [1] 
of superconductivity at 40 K in the previously 
well known material MgB2, which is particularly 
simple from the structural point of view, attracted 
immediately the interest of many research groups 
in the world. Soon after its discovery, several ex- 
periments appeared in the electronic archive cond- 
mat, and subsequently in conventional journals, 
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revealing some microscopic aspects of the super- 
conductivity in this compound that turned out 
later on to be essential for its understanding. The 
boron isotope effect, discovered by Bud'ko et al.[2] 
pointed towards a phonon mediated pairing inter- 
action, and the tunneling spectroscopy, measured 
with a Scanning Tunneling Microscope (STM/S), 
showed a clean BCS density of states with an in- 
triguing low value for the ratio 2A/ksT c = 1.2 
[3], as compared to the value expected from sim- 
ple BCS theory (3.53). First principles investiga- 
tions of the electron-phonon coupling concluded 
that two gaps, related to the two different sets 
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of bands of the Fermi surface would exist in this 
superconductor [4], and that the gap reported in 
Ref.[3] could be associated to one of these sets 
of bands. Within another interesting approach, 
Bascones and Guinea [5] considered the effect of 
the surface as a topological defect in two band 
superconductors to explain the observation of 
a single gap structure in some early tunneling 
experiments[3,6,7,8] . 

Two band superconductivity was theoretically 
proposed by Suhl, Matthias and Walker in 1959[9], 
and has been discussed since then in relation to 
other materials. For instance, important tunneling 
experiments were made by Binnig et al. in 1980 
[10] on Nb doped SrTi03. These authors were able 
to measure in a very neat way the two quasipar- 
ticle peaks in the superconducting tunneling den- 
sity of states due to the opening of two gaps cor- 
responding to the two bands of the Fermi surface 
of this compound, filled one after the other by the 
Nb doping. The results were compared to reduced 
SrTiC>3, where only a broadened single peak struc- 
ture is observed, together with a four times larger 
residual resistivity. Interband scattering due to an 
increasing amount of defects or impurities is indeed 
expected to gradually merge multigap features in 
the density of states into the typical BCS single 
gap curve[ll,12,13]. 

Even if first measurements of the tunneling den- 
sity of states reported the observation of a sin- 
gle gap structure corresponding to different val- 
ues of the superconducting gap (with 2A/ksT c be- 
tween 3 and 4 [6,7,8]), more recent experiments 
[14,15,16,17,18] are in good agreement with the 
model proposed in [4] . Therefore, we can now state 
that tunneling spectroscopy experiments made in 
this material support well the two gap scenario. As 
it is well known, other very powerful macroscopic 
thermal and magnetic measurements, as well as 
microscopic spectroscopic probes, which will be 
treated in detail by other authors in this volume, 
produced more and more compelling evidence in 
this direction (see e.g. [19,20,21,22,23,24,25,26]). 

Previously existing band structure calculations 
[27] were reconsidered in order to understand 



the microscopic nature of the superconducting 
state [28, 29, 30]. Precise calculations based on the 
strong coupling formalism of Eliashberg have 
been recently made by the authors of Ref.[31]. 
It has been shown, and confirmed by comparing 
de Haas van Alphen experiments with the band 
structure calculations [32,33], that the electron 
phonon interaction varies strongly on the Fermi 
surface. Mainly four sheets are found, two of them 
being near cylinders due to the two dimensional 
antibonding states of boron p xy orbitals, and two 
three dimensional sheets from the it bonding and 
antibonding states of the boron p z orbitals[32,34]. 
The superconducting gap A( k ) has been calcu- 
lated as a function of temperature and magnetic 
field and is non zero everywhere on the Fermi sur- 
face. However, its size changes strongly clustering 
around two main values related to the a and 7r 
sheets[31]. The smaller value, which was measured 
in our first experiments [3], appears to be related 
to the three dimensional n sheets, and the larger 
one to the two dimensional a sheets. 

We can therefore affirm that today it is widely 
accepted that MgB 2 is a multiband superconduc- 
tor where the basic physics can be theoretically un- 
derstood in a very sound way. This opens a new and 
interesting challenge to improve the accuracy of the 
experimental results and push the theoretical de- 
velopments towards more and more narrow limits. 
The search for new phenomena associated with two 
band superconductivity appears especially promis- 
ing. With respect to other known multiband super- 
conductors, as e.g. doped SrTiO3[10] or the nickel 
borocarbides[35], MgB 2 has the considerable ad- 
vantage of having a high critical temperature and 
a relatively simple Fermi surface. 

Local vacuum tunneling spectroscopy at low 
temperatures, made possible with the use of 
the STM/S, is a very promising technique, as it 
gives a direct measurement of the local density 
of states [36]. This is even more useful in the case 
of materials, as MgB 2 , where the superconduct- 
ing gap is strongly k dependent. The use of an 
atomically sharp counter-electrode whose (x,y,z) 
position can be changed at will in-situ, gives in 
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principle access to detailed information about the 
superconducting density of states in very differ- 
ent regions of the Fermi surface. However, to get 
fully reliable data we need to reduce as much as 
possible all pair breaking effects acting on the 
surface of the condensate. Careful filtering of all 
wires arriving to the sample turns out to be ex- 
tremely important. The measured spectra most 
accurately approach the local superconducting 
density of states (LDOS) of the sample when 
the measurement is made near T=0K to reduce 
thermal smearing. Obviously, characterizing and 
controlling the surface of the sample is of extreme 
importance for a proper use of this technique. 

We should remark that, although other type 
of spectroscopic measurements provide very use- 
ful information about the density of states, the 
STM/S has the advantage of being a direct probe 
where a very high resolution in energy can be 
achieved by measuring at sufficiently low temper- 
atures. For example, photoemission spectroscopy 
also gives a direct measurement of the density of 
states in different parts of the Fermi surface, but 
present experiments [26] have much less resolution 
in energy than STM/S. Point contact spectroscopy 
has been widely used in the case of MgB 2 (see e.g. 
[19,20]), and it is also a very powerful, but less 
direct probe of the density of states. 

Recently, small single crystals of MgB 2 have 
been grown in a few laboratories[37,38,39], and 
some results on tunneling spectroscopy have been 
reported[18]. At zero magnetic field, STS in the c 
direction at low temperatures has been found in 
agreement with our early result on single grains 
of this material[3]. However, to get a proper fit to 
BCS theory a pair breaking parameter and a finite 
gap distribution, which both account for experi- 
mental smearing, need to be introduccd[18]. The 
observation of a hexagonal vortex lattice is a new 
result, although more experiments are needed to 
understand the microscopic aspects of the Shub- 
nikov phase in this multiband superconductor[18]. 
Other STM/S experiments at 4.2 K in oriented 
thin films, detailed in another contribution to this 
volume, have also found two gap features associ- 



ated to different tunneling directions [6]. 

Several questions remain open which need to be 
addressed using low temperature STM/S. For ex- 
ample, given that this material is a multiband su- 
perconductor, how do the multiple gaps reflect in 
the tunneling characteristics at a given point on 
the surface? How can we extract information about 
the dominant gap structure in this material? How 
shall we deal with and identify the effect of inter- 
band scattering? For instance, would it be possible 
to find a region of the surface where a single gap is 
observed? A zero bias peak due to Andreev bound 
states, which could be related to complex phase re- 
lations between the Cooper pair wavefunction in 
the different bands, has been predicted in multi- 
band, s-wave, superconductors [40]. Are there hints 
of such a peak in STS experiments? On the other 
hand, S-S junctions (which can sometimes be found 
in STM/S experiments) in conventional single gap 
superconductors show a series of small features in 
the conductance below the superconducting gap 
due to multiple Andreev reflection processes (sub- 
harmonic gap structure). At which voltage should 
we expect these features to appear in two gap su- 
perconductors? Actually, how do multiple Andreev 
reflections occur between two multiband supercon- 
ductors?. 

In this paper we present some of our STM/S re- 
search in MgB 2 , made after the first results pub- 
lished in [3] , and try to give insight into some of the 
above mentioned opened questions. New data show 
in both, single crystals and small grains, the same 
reproducible behavior independent of the sample 
preparation method. The two gap structures are 
identified in most typically found spectra. We also 
discuss some serendipitous events where a piece of 
MgB 2 accidentally attached to the tip permits the 
observation of clear S-S tunneling features, includ- 
ing the first (to the best of our knowledge) obser- 
vation of the subharmonic gap structure. 
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2. Experimental and sample preparation 

All cryostats used to measure the data presented 
here (a 4 He, a 3 He and a dilution refrigerator) are 
equipped with similar, home built, STM set-ups 
and electronics which have been carefully filtered 
to produce a clean electromagnetic environment in 
the STM chamber. Measurements in Pb, Al and 
NbSe2 were used to test the energy resolution of all 
set-ups[41,42]. All STM's feature a very high me- 
chanical stability and macroscopic in-situ position- 
ing capabilities (x-y table) in a 2x2 mm 2 window. 
We have measured single crystals and small grains 
prepared in different ways. The single crystals were 
grown under pressure as described in [37] and are 
needle like with typical sizes of 0.1x0.06x0.7 mm 3 
and the c axis along its long direction. Samples 
were introduced (following its long axis) in a small 
hole made on a sample holder (a flat gold surface) 
and glued with a small drop of conductive silver 
epoxy, which was previously put into the small 
hole. About half of the samples remained outside 
the hole. Once a sample was mechanically fixed to 
the sample holder, a small clean blade was used to 
cut it, and the tip was positioned (with the help of 
an optical microscope and the in-situ positioning 
capabilities of our STM) on top of the remaining 
freshly broken surface. The whole set-up was then 
cooled as fast as possible. The tunneling conditions 
were always good, leading to reproducible images 
and spectra over the whole surface. However, the 
topography is complex, showing a large amount 
of inclined small terraces at the nanoscopic scale, 
observed with STM and AFM, indicating that we 
could not find any cleaving plane. SEM images at 
a larger scale also showed a terraced surface. This 
means that the orientation of the surface at the 
length scale of the STM experiment is unknown 
and does not coincide with the breaking plane. The 
small grains were prepared following the procedure 
described in Rcf. [3]. After dispersing them in an 
acetone bath in ultrasound, a small drop was de- 
posited onto a flat gold surface, the liquid was evap- 
orated at 80°C and the grains were pressed into the 



gold surface using a small synthetic ruby. A subse- 
quent bath in acetone and ultrasound removed the 
loosely fixed grains. Following this process, some 
of the samples were heated in a Ar atmosphere of 
600mbar at 200° C during two hours to promote 
gold diffusion and improve the mechanical fixing 
of the grains. 

We made an extensive topographical characteri- 
zation of the surface at all locations where we took 
conductance versus bias voltage curves. However, 
here we report only about the most characteristic 
measured spectra (STS) , because we cannot clearly 
associate the features of the spectra presented here 
to specific features of the surface topography. 



3. Results and discussion 

3.1. 7r and a bands 

In this paper we report about many experi- 
ments made in macroscopically different places 
of a large number of samples, prepared in differ- 
ent conditions. The goal is to show the behavior 
which, to our understanding, best reflects the in- 
trinsic properties of this superconductor. Millions 
of tunneling conductance versus bias curves were 
measured. Most of them were taken in points with 
good tunneling conditions (reproducible imaging 
as a function of the bias voltage and work func- 
tion of several tenths of eV to several eV), and 
the results (measured in all samples and prepara- 
tion conditions) are comprised within two limiting 
behaviors. Before discussing them we should re- 
mark that an STM/S experiment always probes 
a part of the Fermi surface, depending on the 
particular atomic configuration of tip and sample 
surface [36, 43]. Therefore, a multiband material as 
MgB 2 is expected to show a wide variety of be- 
haviors, because, depending on the location, one 
will probe different parts of the Fermi surface. 

In Fig.l we show the tunneling conductance ob- 
tained in a small grain and in a single crystal re- 
spectively. These curves show clearly the partial 
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Fig. 2. Tunneling spectroscopy curves (circles) of MgB2 
obtained in a single crystal (a) and in a grain (b) 
(Rat = 1MQ), together with the distribution of values of 
the superconducting gap used to fit the experiment (lines) . 
At this precise location, the electrons contributing to the 
tunneling current come both from the 7r and a bands. 

is of no use to obtain an adecuate fit to the exper- 
iments. 

To obtain an idea of the contribution to the 
tunneling current from the different bands, we 
phcnomcnologically model the measured super- 
conducting density of states by the sum over a 
normalized distribution of BCS densities of states 
corresponding to a distribution of values of the 
superconducting gap A. As shown in Fig.l, even 
when we deal with the simplest and most fre- 
quent situation, i.e. a single gap structure around 
1.9mV corresponding to the 7r band, a single val- 
ued gap is not enough to obtain a the best fit to 
the experiment. The observed quasiparticle peaks 
are not as high as expected from most simple s- 
wave BCS theory. They are smaller because of a 



Fig. 3. Different spectra taken at different positions on a 
single grain (circles), and its associated fits (lines and in- 
set). In (a) we observe a single, broadened quasiparticlc 
peak, whereas (b) and (c) show intermediate situations be- 
tween those presented in figs.l and 2. We associate this 
behavior, which is also observed in single crystalline sam- 
ples, to interband scattering effects at the surface of MgB2 
(see text). 

increasing the smaller gap, decreasing the larger 
gap and leading therefore to a smaller critical 
temperature[ll,12,13,30]. However, even bad qual- 
ity samples of MgB2 do not show big changes in 
the superconducting critical temperature, and also 
show both superconducting gaps with values which 
do not differ significantly from the best samples. 
The theoretical work in Rcf.[13] deals with this ap- 
parent contradiction by proposing that interband 
scattering is particularly small in this compound. 
It would be highly desirable to obtain experimen- 
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Fig. 4. Subharmonic gap structure observed in spectra char- 
acteristic for tunneling between two single grains of MgB2, 
taken at 2K. The inset shows a zoom over the relevant 
voltage region with arrows at each peak corresponding to a 
submultiple of 2A n , and at the peak originating in the con- 
ductance due the Joscphson effect (arrow to the left). Note 
that the measured negative resistance around zero bias is 
a characteristic feature of junctions where the Joscphson 
effect is observed. 

are still observed. 

3.3. S-S junctions 

As described earlier, in the experiments made in 
single grains, we use the x-y table to change macro- 
scopically the position of the tip with respect to the 
sample. This enables us to clean the tip on the gold 
surface and then to measure STM/S with a normal, 
clean gold tip on a single MgB2 grain. Sometimes 
however clear S-S tunneling spectra appeared, due 
to a loosely connected single grain which was at- 
tached by accident to the tip apex, as already men- 
tioned in [3], where we also reported that the ob- 
served behavior was in agreement with the value 
of A found with a clean normal tip. 

In Fig. 4 we represent conductance curve ob- 
tained in one of those experiments [3]. A peak in 
the conductance due to the Joscphson current 
neatly appears at zero bias, and the so-called 
subharmonic gap structure, characteristic of S-S 
junctions (inset). It is indeed well known that mul- 
tiple Andreev reflections occurring at the junction 
provide a vehicle for a finite conduction within the 



Fig. 5. The temperature dependence of the S-S spectra as 
measured in a junction made of two parts of single crys- 
talline MgB2 (curves displaced by 1 in the y axis). Dashed 
lines are fits to these curves. For the lowest temperature 
curve we use a distribution of values of the superconduct- 
ing gap as shown in the inset at the bottom of the fig- 
ure. For higher temperatures, this distribution has to be 
changed slightly due to the temperature variation of the 
superconducting gaps. The contribution from a band elec- 
trons to the tunneling current is small and leads to a broad 
peak, which has been zoomed out in the figure of the inset. 
Straight lines at the bottom of each curve indicate zero 
conductance for each temperature. 

of taking the derivative of the conductance 
function of the bias voltage d(^)/dV [49], which 
enhances the features due to the dominating gaps 
in the distribution, we can best identify the char- 
acteristic energy corresponding to the physical 
process which gives these peaks, as a maximum 
(minimum) in d(^)/dV for positive (negative) 
voltages. The features at 3.8 and 9.5 mV are 
clearly due to conventional tunneling processes 
and are also obtained from the fit to the experi- 
mental curves shown in Figs. 4 and 5. The feature 
in Fig. 6a at 3.8mV (corresponding to the highest 
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Fig. 6. The derivative of the conductance as a function of 
the bias voltage in the S-S junction at 2K(a). This type of 
plots is used to identify the voltage position of the main 
peaks obtained in the conductance (Figs. 4 and 5). Note 
that the voltage range is much larger than the one in inset 
of Fig. 4, so that the peaks of the subharmonic gap structure 
of 2A W corresponding to n > 2 (located at 2A n /n, see 
right bottom inset of Fig. 4) appear as a very small feature 
at the center of this plot. The arrows show the peaks which 
are best identified in this voltage range. In (b) we show 
the temperature dependence of A CT and A^ as determined 
by using plots as (a). 

peak in Figs. 4 and 5), is located at two times the 
dominant gap in the distribution corresponding 
to the 7r band gaps, 2A„. The contribution to the 
tunnelling current from the a band electrons is too 
small in these junctions to give a feature at 2A CT 
as shown by the fit in Fig. 5 and its associated gap 
distribution. The clear feature observed at 9.5mV 
is located at the sum of the dominant gaps of each 
set of bands, A n + A a . The physical origin of this 
effect is that, when the bias voltage between both 
electrodes coincides with the sum of the peaks of 
both distributions A T + A CT , the high peak in the 
density of states corresponding to the tt band gaps 
enhances the small bump corresponding to the a 



band gaps and gives a clearly measurable peak in 
the conductance. 

From Fig. 6a we can also obtain the location of 
the peaks related to the subharmonic gap structure 
associated to the A T gap, already discussed in con- 
nection with Fig. 4 (in Fig. 6a, for clarity, we only 
highlight 2A 7T /n with n = 2). Interestingly, we can 
also observe another feature at A a , possibly related 
to multiple Andreev reflections which enhance the 
conductance at the voltage corresponding to the 
n = 2 subharmonic gap structure associated to the 
a band gaps. This point should be analyzed in more 
detail by extending presently available theoretical 
models to multiband materials. To the best of our 
knowledge, no other previous measurement of this 
characteristic feature of S-S tunnel junctions has 
been made in multigap superconductors. 

By using the same plot for higher temperature 
curves, we can also easily obtain the temperature 
dependence of the dominant superconducting gap 
features at the tt and a bands. The result is repre- 
sented in Fig. 6b. Note that, while A w stays fairly 
constant, A CT changes more rapidly as a function 
of temperature. 



4. Conclusion 

We have presented new results on Scanning Tun- 
neling Spectroscopy in single crystals and grains 
of MgB2- The tunneling spectra clearly reflect the 
multiband character of superconductivity in this 
material, and are well fitted by two band BCS ex- 
pressions. We mostly find features associated to the 
7r band gap, which clearly dominates the contri- 
bution to the tunneling spectra, although features 
due to the a band gaps were also detected. Both 
bands have a non negligible distribution of values 
of the superconducting gap, in agreement with the 
recent theoretical calculations of [31]. 

In some locations of the surface we can observe 
interesting effects with which we can extend our 
present understanding of this material. Interband 
scattering effects lead to the observation of a sin- 
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gle gap structure. S-S junctions show clearly the 
Josephson effect and the appearance of peaks at 
sub-multiples of the superconducting gap charac- 
teristic of this type of junctions (subharmonic gap 
structure). 

Clearly, many questions remain open which need 
to be addressed in future experiments. The obser- 
vation of atomic resolution images of the surface 
is necessary to understand the dominant role of 
the 7r band gap in our measurements. It is indeed 
striking that other important experiments, as e.g. 
NMR or Raman scattering [25,51], preferentially 
observe features associated to the a band gap. In 
NMR [51], even no trace of the 7r band gap has 
been measured. Although these last experiments 
were made at a relatively high magnetic field (IT), 
which could have influenced in some way the ir 
band gap, the reasons for the observation of one 
or the other gap structure in different experiments 
needs to be addressed. Furthermore, the search for 
phonon features in the electronic tunneling density 
of states in order to obtain more information about 
the pairing interaction is an open issue. It should 
also be interesting to study in more detail S-S junc- 
tions of this material in a situation with a larger 
contribution from the a band electrons, as for in- 
stance the one presented in fig. 2b, to the tunneling 
current to search for new features in the subhar- 
monic gap structure and in the Josephson effect. 
The resulting curves should be highly non linear, 
opening the way to new interesting applications. 

Theoretical calculations pointing out possible 
new phenomena to be observed in N-S and S-S tun- 
neling experiments associated to the multigap na- 
ture of superconductivity in this compound would 
be very stimulating. In any case, MgB 2 is clearly a 
very intriguing superconductor with a remarkably 
simple Fermi surface and high critical temperature. 
This has made, and will continue to make, feedback 
between theory and experiment especially fruit- 
ful. We can therefore expect new and surprising 
physics to emerge from this interaction in the fol- 
lowing years. 



Acknowledgement We specially acknowledge F. 
Guinea and E. Bascones for very fruitful discus- 
sions. We also acknowledge J. Ramirez, N. Agrait 
and J.L. Martinez. Support from the ESF pro- 
gramme Vortex Matter in Superconductors at Ex- 
treme Scales and Conditions (VORTEX), as well as 
from MCyT (MAT2001-1281-c02-01, Spain), the 
Comunidad Autonoma de Madrid (Spain) and the 
European Social Fund is acknowledged. The sin- 
gle crystals of MgB 2 were prepared in support by 
the New Energy and Industrial Technology De- 
velopment Organization (NEDO) as Collaborative 
Research and Development of Fundamental Tech- 
nologies for Superconductivity Applications. 



References 

[I] J. Nagamatsu, N. Nakagawa, T. Muranaka, Y. Zenitani, 
J. Akimitsu, Nature (London), 410, 63 (2001). 

[2] S.L. Bud'ko, G. Lapertot, C. Petrovic, C.E. 
Cunningham, N. Anderson, P.C. Canficld, Phys. Rev. 
Lett. 86, 1877 (2001). 

[3] G. Rubio-Bollinger, H. Sudcrow, S. Vieira, Phys. Rev. 
Lett. 86, 5582 (2001). 

[4] A. Y. Liu, I. I. Mazin, and J. Kortus Phys. Rev. Lett. 
87, 087005 (2001) 

[5] E. Bascones and F. Guinea, Phys. Rev. B 65, 174505 
(2002) 

[6] G. Karapctrov, M. Iavarone, W. K. Kwok, G. W. 
Crabtrce, and D. G. Hinks Phys. Rev. Lett. 86, 4374- 
4377 (2001) 

[7] H. Schmidt, J. F. Zasadzinski, K. E. Gray, and D. G. 
Hinks Phys. Rev. B 63, 220504 (2001) 

[8] A. Sharoni, I. Felncr, and O. Millo, Phys. Rev. B 63, 
220508(R) (2001) 

[9] H. Suhl, B.T. Matthias and L.R. Walker, Phys. Rev. 
Lett., 3, 552 (1959) 

[10] G. Binnig, A. Baratoff, H. E. Hocnig, and J. G. 
Bednorz Phys. Rev. Lett. 45, 1352-1355 (1980) 

[II] C.C. Sung, V.K. Wong, J. Phys. Chem. Solids, 28, 
1933 (1967). 

[12] A. A. Golubov and I. I. Mazin, Phys. Rev. B 55, 15146- 
15152 (1997) 



10 



[is; 
[i4; 

[15 



I. I. Mazin, O. K. Andersen, O. Jepsen, O. V. Dolgov, 
J. Kort us, A. A. Golubov. A. B. Kuz'menko, D. van der 
Marel, |:ond-mat/0204013 



H. Schmidt, J. F. Zasadzinski, K. E. Gray, and D. G. 
Hinks Phys. Rev. Lett. 88, 127002 (2002) 

M. Iavarone, G. Karapetrov, A. E. Koshele y. W 
K. Kwok. G . W. Crabtree, D. G. Hinks, |cond- 



mat 



[ie; 

[i7; 
[is; 

[19 

[20 

[21 
[22 

[23 
[24 
[25 



F. Giubileo, D. Roditchev, W. Sacks, R. Lamy, D. X. 
Thanh, J. Klein, S. Miraglia, D. Fruchart, J. Marcus, 
and Ph. Monod Phys. Rev. Lett. 87, 177008 (2001) 

M. H. Badr, M. Freamat, Yu. Sushko, and K.-W. Ng 
Phys. Rev. B 65, 184516 (2002) 

M. R. Eskildsen, M. Kugler, S. T anaka. J. Jun. S. M . 
Kazakov, J. Karpinski, O . Fischer, pond-mat /0207394 

P. Szabo, P. Samuely, J. Kacmarck, T. Klein, J. 
Marcus, D. Fruchart, S. Miraglia, C. Marcenat, and A. 
G. M. Jansen Phys. Rev. Lett. 87, 137005 (2001) 

Yu.G. Naidyuk, I.K. Yanson, L.V. Tyutrina, N.L. 
Bobrov, P.V. Chubov, W. N. Kang. Hveong-J in Kim, 
Eun-Mi Choi, Sung-Ik Lee, 



/0203329| . 



:ond-mat/0112452 



F. Bouquet, R. A. Fisher, N. E. Phillips, D. G. Hinks, 
and J. D. Jorgensen Phys. Rev. Lett. 87, 047001 (2001) 

H. D. Yang, J.-Y. Lin, H. H. Li, F. H. Hsu, C. J. Liu, 
S.-C. Li, R.-C. Yu, and C.-Q. Jin, Rev. Lett. 87, 167003 
(2001) 

Y. Wang, T. Plackowski, A. Junod, Physica C, 355, 
179 (2001). 

F. Manzano, A. Carrington, N.E. Husscy, S. Lee and 
A. Yamamoto, Phys. Rev. Lett. 88, 047002 (2002) 

J. W. Quilty, S. Lee, A. Yamamoto, and S. T aiima. 
Phvs. Rev. , Lett. 88, 087001 (2002) and |cond- 



mat/0206506 



[26] S. Tsuda, T. Yokoya, T. Kiss, Y. Takano, K. Togano, 
H. Kito, H. Ihara, and S. Shin Phys. Rev. Lett. 87, 
177006 (2001) 

[27] D.R. Armstrong and P.G. Perkins, J. Chem. Soc. 
Faraday Trans. 275, 12 (1979). 

[28] J. Kortus, I.I. Mazin, K.D. Belashenko, V.P. Antropov 
and L.L. Boyer, Phys. Rev. Lett 86, 4656 (2001). 

[29] J.M. An, W.E. Pickett, Phys. Rev. Lett 86, 4366 
(2001). 

[30] A. Brinkman, A. A. Golubov, H. Rogalla, O. V. 
Dolgov, J. Kortus, Y. Kong, O. Jepsen, and O. K. 
Andersen Phys. Rev. B 65, 180517 (2002) 

[31] H.J. Choi, D. Roundry, H. Sun, M. L. Cohen, S.G. 
Louie, Nature, 418, 758 (2002) 

[32] E. A. Yelland, J. R. Cooper, A. Carrington, N. E. 
Hussey, P. J. Meeson, S. Lee, A. Yamamoto and S. 
Tajima, Phys. Rev. Lett. 88, 217002 (2002) 



[33 

[34; 

[35 
[36 

[37 
[38 

[39 

[10 
[41 

[42 
[43 
[44 
[45 
[46 

w. 

[48 

[49 
[50 
[51 



I. I. Mazin and J. Kortus Phys. Rev. B 65, 180510 
(2002) 

H. Uchiyama, K.M. Shen, S. Lee, A. Damascelli, D.H. 
Lu, D.L. Fheng, Z.X. Shen and S. Tajima, Phys. Rev. 
Lett. 88, 157002 (2002) 

P.C. Canfield, P.L. Gammcl, D.J. Bishop, Phys. Today 
(1998), 51, pp. 40-46. 

C.J. Chen, Introduction to Scanning Tunneling 
Microscopy, Oxford Series in optical and imaging 
sciences, Oxford University Press (1993). 

S.Lee, T.Masui. H.Mori. Yu.Eltsev. A. Yamamoto, 



S.Tajima, xmd-mat/0207247 



J. Karpinski, M.Angst, 
J. Jun, S.M. Kazakov, R.Puzniak, A.Wisnicwski, J.Roos, 
H.Keller, A. Perucchi, L. De giorgi. M. Eskildsen . 
P.Bordet, L.Vinnikov, A.Mironov, cond-mat/0207263| . 

| Y. Machida, S. Sasa ki. H. Fuiii. M. F uruvama. I. 
Kakeya, K. Kadowaki, |cond-mat /0207658[ 



K. Voelker, M. Sigrist, cond-mat/0208237 



H. Suderow, M. Crespo, P. Martinez-Samper, J.G. 
Rodrigo, G. Rubio-Bollinger, S. Vieira, N. Luchier, J. P. 
Brison, P.C. Canfield Physica C, 369, 106 (2002). 

H. Suderow, E. Bascones, A. Izquierdo, F. Guinea, S. 
Vieira, Phys. Rev. B, 65, 100519 (R) (2002). 

V.M. Silkin, E.V. Chulkov and P.M. Echenique, Phys. 
Rev. B, 64, 172512 (2001). 

R.C. Dynes, V. Narayanamurti, J. P. Garno, Phys. Rev. 
Lett., 41,p. 1509 (1978). 

H.F. Hess, R.B. Robinson, J.V. Waszcak, Phys. Rev. 
Lett. 64, p. 2711 (1990). 

H. Suderow, P. Martinez-Samper, N. Luchier, J. P. 
Brison, S. Vieira and P.C. Canfield, Phys. Rev. B, 64, 
R020503 (2001). 

P. Martinez-Samper, H. Suderow, S. Vieira, N. 
Luchier, J. P. Brison, P. Lejay, S. Vieira, submitted. 

Y. Wang, F. Bouquet, I. Sheikin, P. Toulemonde, B. 
Rcvaz. M. Eistcrer, H. W. Weber, J. Hinderer, A. Junod, 



:ond-mat/020816£ 



E.L. Wolf, "Principles of Electron Tunneling 
Spectroscopy", Oxford University Press (1989). 

Y. Zhang, D. Kynion, J. Chen, J. Clarke, D.G. Hinks, 

G. W. Crabtree, Appl. Phys. Lett., 79, 3995 (2001). 

H. Kotegawa, K. Ishida, Y. Kitaoka, T. Muranaka, and 
J. Akimitsu, Phys. Rev. Lett. 87, 127001 (2001) and 

H. Kotegawa, K. Ishida, Y. Kitaoka, T. Muranaka, N. 
Nakagawa, H. Takagiwa and J. Akimitsu, Phys. Rev. B 
66, 064516 (2002) 



11 



